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1 Introduction 



In recent years there has been huge progress in exactly solving the A/" = 4 
Super- Yang-Mills theory in the planar limit PP-[IS] • This is particularly 
remarkable as = 4 SYM is an interacting four-dimensional gauge theory 
with highly nontrivial dynamics. All previous examples of exact solvability in 
Quantum Field Theories were restricted either to two dimensional theories or, 
in the case of higher dimensional theories, to some supersymmetric subsector 
or to theories with much simpler dynamics like topological field theories. 

The solvability of A/" = 4 SYM arises due to the AdS/CFT correspondence 
fiQ[ [T71 [18] , according to which the gauge theory is equivalent to string theory 
in AdS^ X 5*^. Therefore Af = 4 SYM avoids the no-go theorems for inte- 
grable field theories in more than two-dimensions by translating its dynamics 
into properties of the two-dimensional worldsheet QFT of the superstring in 
AdS^ X which is an integrable QFT. 

Since N' = 4 SYM is a conformal field theory, all correlation functions 
of local operators are, in principle, determined by a much smaller set of 
data: the set of conformal dimensions (equivalently anomalous dimensions 
of gauge theory operators) and the OPE coefficients. These data can be 
extracted from the knowledge of 2- and 3-point correlation functions. 

In the AdS/CFT case, however, the anomalous dimensions are extracted 
not from 2-point correlation functions but directly as eigenvalues of the di- 
latation operator, which translates to energies of string states in AdS^ x S^. 
Therefore the problem of finding all anomalous dimensions reduces to find- 
ing the energy levels of the two-dimensional worldsheet QFT on a cylin- 
der. Currently we have a very complete understanding of the spectrum of 
conformal dimensions which is described by a set of Thermodynamic Bethe 
Ansatz equations [II]-[I5]. The methods used are similar to the ones em- 
ployed for relativistic integrable field theories, although their generalization 
to the AdS / CFT case is far from trivial due to many unique features of the 
worldsheet QFT. 

For the case of OPE coefficients, however, there seems to be no alternative 
to a direct computation of 3-point correlation functions. It is convenient to 
classify the operators into three main groups, depending on their behaviour 
at strong coupling: 'light', 'medium' and 'heavy' operators. The 'light' oper- 
ators are BPS and dual to supergravity fields. Consequently their anomalous 
dimensions do not depend on the coupling. The next class of operators are 
the lightest massive string states whose dimensions scale like Aj. A classical 
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example of these 'medium' operators is the Konishi operator. Finally, the 
'heavy' operators have large charges (of the order of A2) and are dual to 
classical string states with anomalous dimensions scaling like As [TU]-|21]. 
Although very useful, this is only a rough and nonexhaustive classification. 
BPS operators with large charges may for all practical purposes behave like 
'heavy' operators. There may be operators with dimensions like A 2 which 
may be very quantum and without a classical string description (like BPS 
operators with two large charges). 

The techniques for computing 3-point correlation functions are very well 
developed for the case of 'hght' operators, i.e. BPS operators dual to super- 
gravity fields [22]-[23]. Unfortunately these OPE coefficients are protected 
and do not depend on gauge theory coupling. 

For unprotected operators, the techniques for computing even 2-point cor- 
relation functions have been developed only recently [24j (but see also [22] 
and [26j). These results have been used to compute OPE coefficients between 
two 'heavy' and one 'light' operator using the known classical solution cor- 
responding to a 2-point function and integrating a supergravity propagator 
over the classical string worldsheet [27], [28], [29]-[l0]. It has been further 
extended for correlators involving two Wilson loops and a 'light' operator 

mm 

An intermediate case recently considered in the literature involved a 
geodesic approximation for the three operators [43J. Such an approxima- 
tion may be very relevant to the case of three 'medium' operators which are 
not sufficiently heavy to generate an extended, non-pointlike, surface. 

The goal of this paper is to compute (the AdS part of) the 3-point corre- 
lation functions of three 'heavy' operators. We assume that these operators 
do not have any spin in AdS^. The main difficulty lies in the fact that a 
novel type of a classical solution has to be found. Moreover, in contrast to 
the spectral problem, there is no analog of this problem for conventional rel- 
ativistic integrable field theories, therefore we do not have any guide in this 
respect. 

The computation of the OPE coefficients for three 'heavy' operators is 
especially interesting in view of the fundamental importance that the in- 
tegrable classification of finite-gap spinning string solutions and their com- 
parision with 1-loop Bethe Ansatz results had in arriving at the all-loop 
interpolation. We hope that a similar comparision with weak coupling data 
[1^ H5] will be very illuminating also in the case of OPE coefficients. 
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The plan of this paper is as follows. In section 2, we briefly review the 
case of 2-point correlation functions, and in section 3 the general features of 
the problem of finding 3-point correlation functions. In section 4, we give 
an overview of our approach to this problem, in order for the reader to not 
get lost in the technicalities. In section 5, we review Pohlmeyer reduction, 
and give our main technical results necessary for later computations — we 
solve functional equations for the products between the solutions of the linear 
system on the 3-punctured sphere and give formulas for reconstructing the 
classical solution in AdS2 from the Pohlmeyer data. We then proceed, using 
these results, to evaluate in the next 3 sections the two main parts of the AdS 
contribution to the correlation function. We give our final result in section 9 
and discuss the limits of small and large anomalous dimensions and the link 
of the latter to the Painleve III transcendent. We close the paper with a 
discussion and several appendices with some technical details. 

2 Two-point correlation functions 

In this section we will briefly review the computation of 2-point correlation 
functions for the class of operators that we are considering in this paper, 
namely operators corresponding to classical string solutions with no charges 
in AdS |I9]-[2T]. 

The approach introduced in [21], amounts to computing a cylinder am- 
plitude (with additional wavefunctions included in order to project on the 
specific string state that we are interested in) with the boundary conditions 
such that the string worldsheet approaches two given points on the bound- 
ary regularized by a cut-off z = S. This is done by a classical computation, 
where the corresponding solution is just a geodesic in the AdS part and in 
the S^ part coincides with the unmodified S^ spinning string solution used in 
the conventional calculation of the anomalous dimension. Then one performs 
a saddle point evaluation of the integral over the modular parameter. The 
outcome i^ that the saddle point value of the modular parameter is purely 
imaginary thus effectively making the worldsheet Euclidean. 

This generic pattern indicates that we could have started directly at the 
saddle point, with the worldsheet having the topology of a 2-punctured sphere 
(again with small disks corresponding to z = S cut out), and the Euclidean 
solution satisfying Virasoro constraints. 

^See specific examples in |24j. 



4 



The Euclidean solution for the operators in question has the following 
simple structure. The AdS^ part reduces just to a geodesic in the AdS2 
subspace which contains the two gauge theory operator insertion points on 
the boundary. Explicitly we have 

R R 1 /n -1 \ 

X — — tanh KT + xo z — ; (2.1) 

2 2 cosh/tr ^ ^ 

where the distance between the operator insertion points is Xi — X2 = R. 
Imposing the target space cut-off z = S translates into a worldsheet cut-off 
which limits the range of r to 

At = - log f (2.2) 

K t 

The part is just the Wick rotated spinning string solution, a simple ex- 
ample being a circular string with equal angular momenta given in terms of 
the standard angular coordinates on C -S"^ by 

ijj = a <Pi = 4>2 — (2.3) 

One should note that due to the i, the solution is inherently complex and in 
fact does not look like any kind of spinning string. Virasoro constraints link 
K and u! through k — vT+u^. In general we have k — A, where A is the 
dimension of the operator. 

The 2-point correlation function is now obtained by i) evaluating the AdS2 
action on the AdS geodesic part and ii) evaluating the action together 
with wavefunction contributions, which transforms the action integral into 
an Euclidean energy integral. 

Exphcitly we get 

exp {-v^^TW log ll -exp {-v^v^lW log |j (2.4) 

^ V ' ^ V ' 

AdS action energy 

which reproduces the 2-point correlation function with the correct value of 
the anomalous dimension 

R 



{0(x,)0{x2)) = ( -J (2.5) 



The wavefunction factors in the part are crucial for the correct answer. 
The Virasoro constraint sets the only free parameter k in the AdS2 geodesic 
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solution (2.1) to be equal to A, the (anomalous) dimension of the corre- 
sponding operatoij^ Therefore, the AdS part of the solution is completely 
determined just by the dimension of the operator in question and not by any 
further details of the specific operator. We will see that the same property 
will hold also for 3-point correlation functions. In the following we will use 
complex coordinates which, in the present case, read 

w = e"+^" w = e^-'" (2.6) 

putting the two punctures at u; = and w = oo. 



Finally, let us note that the target space cut-off £ enters (2.5) essentially 
as the normalization of our operators. This is important as in the case of 
3-point correlation functions we have to retain exactly the same normaliza- 
tion of operators as for 2-point functions in order to extract unambiguosly 
the OPE coefficients. This leads to severe difficulties, such as linking the 
worldsheet cut-offs around the 3 punctures to the single target space cut-off 
z = S. This is highly nontrivial due to the lack of an explicit classical ^4^5' 
solution in the case of 3-point correlation functions. If one would adopt a 
different approachj^of using vertex operators [IHl HTj |26] , then the difficulties 
remain but appear in different places. In the vertex operator approach, one 
computes the worlsheet integral over the whole punctured sphere with the 
vertex operator contributions sitting directly at the punctures. For 2-point 
functions one then just neglects infinities. For 3-point functions it is not 
clear how to control possible finite renormalizations. It would be very inter- 
esting to understand quantitatively the precise dictionary between the two 
approaches. 



3 Three-point correlation functions — general 
features 

In order to compute the 3-point correlation function of heavy operators, we 
have to find a classical solution of string equations of motion in Euclidean 
signature with the topology of a sphere with 3 punctures, with the property 
that the solution close to each puncture associated with a given gauge theory 

^More precisely the dimension A is A = \/\A. 
•^which should of course be in the end equivalent 
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operator Ok looks asymptotically like a solution corresponding to a 2-point 
correlation function of the operator Ok- 

The classical bosonic equations of motion in AdS^ x reduce to two 
independent sets of equations, one on 5'^, the other on AdS^., which are 
coupled together only through the Virasoro constraint 

TAds,{w)+Ts^.{w)=Q (3.1) 

where w is the holomorphic worldsheet coordinate, and TacisA''-^') (I'csp. T55) 
is the classical energy-momentum tensor of the AdS^ (resp. S^) part of the 
(7-model. 

For operators which do not have any spins, the AdS^ part of the problem 
greatly simplifies. Then, without loss of generality, wc put the gauge theory 
operator insertion points of all three operators on a single line. Consequently, 
the AdS^ part of the string solution is contained in an (Euclidean) AdS2 
subspace. The problem, however, does not trivialize as we are not looking 
for a minimal surface but have a prescribed nonzero energy-momentum tensor 

TAdS2{'^)- 

Fortunately, the AdS2 energy-momentum tensor can be explicitly ex- 
pressed just in terms of the anomalous dimensions of the three operators en- 
tering the 3-point correlation function. Prom now on we will denote TAdS2{'^) 
by T{w). 

In order to find the explicit form of T{w), recall that the classical solution 
should approach, at the punctures, 2-point solutions which are explicitly 
known. In particular T{w) for the 2-point solutions is given by 

T2-pointiw) = (3.2) 

Therefore at the punctures T{w) should have at most poles of 2"*^ order with 
the leading coefficients determined by the dimension of the operator in- 
serted at w = Wk- Since T{w) is holomorphic and transforms under inversion 
like a (2,0) tensor 

rw^lr(i) (3.3) 

its form is uniquely determined. Without loss of generality we may put the 
punctures at ty = ±1 and w — 00. Then, for the case of equal conformal 
weights A at w = ±1 and Aoo at ty = 00, T{w) is given by 

n-) = ^Z^, where a'=*-^-l (3.4) 



4 (l-m2)2 A5 



OO 
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In the present paper for simplicity we will predominantly consider the above 

symmetric case. 

To summarize, we thus have to evaluate the action 

exp ^ rSr^ci^^ j " (- V Z ^-^^^^^d^v)j (3.5) 

for a classical solution approaching the operator insertion points Xk at 
w — — 1,1,00 subject to the constraint 

^^^l±^=TM (3.6) 

z 

4 Our strategy 

As described in the previous section, the 3-point correlation functions for 
three heavy operators with no spins in AdS^ factorize into a product of an 
and an AdS2 contribution evaluated for a worldsheet with the topology 
of a 3-punctured sphere. Similarly as for 2-point functions we regularize the 
worldsheet by cutting out small disks of radii Si around the punctures which 
are defined by the condition that on their boundaries 

z^S (4.1) 

where z is the AdS radial coordinate in the Poincare patch (2; = is the AdS 
boundary). £ is the target space cut-off which is taken to be very small. It 
is necessary to ensure that z = S around each puncture in order to have the 
same normalization of operators in 2- and 3-point correlation functions so as 
to unambigously extract the OPE coefficients. 

For the AdS2 part, we have to evaluate the action of the classical solution, 
while for the part wc have to include, in addition, contributions from 
the classical wavefunctions of the external states. Therefore, the 3-point 
correlation function is schematically given by 



y/X r ^Polyakov y/\ r y^Polyakov 

^ V ' 

/e\{£J contribution 

Since both exponents have logarithmic divergences around the punctures, 
it is convenient to subtract and add ^T{w)T{w) regularizing the integrals. 
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This yields 



for the AdS2 part and 



- r -c 

TT J-E\'~AdS2 



'TT 



'TT 



/j, ( 55 contribution- \/TT]-^ J,^ 



'TT 



(4.3) 



(4.4) 



for the part. The first terms in the above expressions are now finite 
and can be integrated over the whole punctured sphere, while the explicit 
dependence on the worldsheet cut-offs appears only in the second integral 
with a known integrand. 



In this paper we will compute the contribution (4.3) together with the 



second term in (4.4), leaving the remaining factor 



f^i contribution- 



'TT 



(4.5) 



for further investigation. 
In order to compute 



1 Polyakov 
■'AdS2 



'TT 



(4.6) 



we will use Pohlmeyer reduction |1S], [13] and adapt the methods of 
to evaluate this expression. Firstly, one transforms the above integral into 
an integral of the wedge product of two closed 1-forms on a double cover 
of S. Secondly, one uses Riemann reciprocity (Riemann bilinear identity) to 
express the integral in terms of products of integrals of the 1-forms on certain 
open cycles. Thirdly, one links the above 1-form integrals to the asymptotics 
in the spectral parameter (^ — )■ 0) of appropriate skew products between 
specific solutions (associated with each puncture) of the Pohlmeyer linear 



system. Thus the evaluation of the integral (4.6) is reduced to the knowledge 



of appropriate skew products as a function of the spectral parameter. 
The remaining integral 



'TT 



(4.7) 



can be evaluated analytically in the small £j limit. The main difficulty lies in 
linking the worldsheet cut-offs {si] to the target space cut-off z = £, without 
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an explicit knowledge of the classical solution. To do that we need formulas 
for reconstructing the classical solution from Pohlmeyer data. Fortunately, 
since the classical solution should approach the known solutions for 2-point 
functions close to the punctures, we can get exphcit formulas relating the 
positions of the gauge theory operator insertion points Xk and the target 
space cut-off S to the worldsheet cut-offs {si} in terms of the skew products 
mentioned above, but this time evaluated at ^ = 1. Using this knowledge, 
the (two copies of the) integral (4.7) yield the standard space-time dependent 
part of the 3-point CFT correlation function, as well as a finite contribution 
to the OPE coefficient expressed in terms of the skew products at ^ = 1. 

The skew products between the specific solutions (of the Pohlmeyer linear 
system) associated to each puncture are therefore a key ingredient in the 
evaluation of the 3-point correlation function. We will often refer to these 
chosen solutions as 'elementary solutions'. 

In the following section we introduce the main features of Pohlmeyer re- 
duction, the elementary solutions associated with each puncture and define 
the skew products. Then we derive and solve functional equations for the 
skew products of the elementary solutions as a function of the spectral param- 
eter ^. Finally we state the reconstruction formulas which link the operator 
insertion points and the target space cut-off to appropriate skew products. 

After this preparatory part, we proceed to evaluate the integral (4.6) in 
section 7 and the divergent contribution (4.7) in section 8. Then we put 
together the obtained formulas into the final AdS contribution to the OPE 
coefficient and analyze the limits of large and small anomalous dimensions 
as well as the extremal limit. 

Before we end this overview, let us remark that the same decomposition 
(4.2) could also be performed for a 2-punctured sphere corresponding to a 
2-point correlation function (of course in this case only two wavefunctions 
would appear). Then it turns out that both the 'nontrivial' parts (4.6) and 

(4.5 ) are identically zero. However it is interesting to note that they vanish for 
quite different reasons. The AdS part (4.6 ) vanishes because it is evaluated on 
a trivial classical solution - a point-like string moving along a geodesic. The 
corresponding Pohlmeyer function is just identically zero and consequently 

(4.6) vanishes. On the S^, however, we deal with arbitrarily complicated 
finite-gap solutions of arbitrary genus, which would have a highly nontrivial 
Pohlmeyer description. Yet, the wavefunction contributions transform the 
classical action into an integral of the energy density (in an appropriate 
coordinate system) and the resulting contribution exactly cancels the 
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integral of vTT. It is tempting to speculate that a similar simplification 
may occur for the case of 3-point functions. 

5 Pohlmeyer reduction 

Contrary to the well known case of Pohlmeyer reduction for minimal surfaces 
in AdS^ [H], [l9], [50], we need to perform Pohlmeyer reduction for classical 
solutions in AdS2 but with a prescribed nonzero energy-momentum tensor. 
Thus the classical solutions in AdS2 are of course not minimal surfaces. On 
the other hand, the full string solution, which takes into account both AdS2 
and 5*^ contributions, is a minimal surface. 

The Pohlmeyer reduction for this case amounts to defining 7(u7, w) through 

^^^^ + ^"^" = cosh 7 (5.1) 

where T is the energy-momentum tensor T{w). Then ^{w,w) satisfies a 
modified form of Sinh-Gordon equation 

dd^=\/Tfsmhj (5.2) 

The solution corresponding to a 2-point function is just ^{w,w) = 0. Con- 
sequently, the boundary conditions close to each puncture are 

7^0 (5.3) 

For the case relevant to 3-point correlation functions, T{w) has two zeroes, 



and thus the form of Pohlmeyer reduction given by (5.1) is inconvenient as 
it would imply that all first derivatives vanish 

dz = dz = dx = dx = (5.4) 

at the zeroes of T{w). This would be a very nongeneric situation as each 
such single equation gives a codimension one subspace. Their intersection 
is generically empty. This is even the case for pointlike strings appearing in 
2-point functions. Consequently we will assume, as is the case for polygonal 



Wilson loops, that the right hand side of (5.1) is everywhere nonzero. This 



implies that 7 has to have logarithmic singularities at the zeros of T{w) 
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To avoid this drawback, it is convenient to redefine 7 through 

7 = 27 log TT (5.5) 



Now (5.1) takes the form 

dxdx + dzdz 



(e27 ^ TTe-^^) (5.6) 



which does not lead to any problem at the zeroes of T[w). The equation of 
motion becomes 

dd-i = ^ (e'^ - (5.7) 



This is virtually the same as the setup for Wilson loop [50], [51] but with 
the polynomial defining the polygonal Wilson loop substituted by T{w). We 
will discuss the similarities and differences in more detail at the end of the 
present section. 

It is well known that the modified sinh-Gordon model is integrable. It is 
easiest to verify by making a holomorphic change of worldsheet coordinates 
to map this model into ordinary sinh-Gordon. However, due to the rather 
complicated analytical structure of the resulting domain we will not use this 
mapping in the sequel. 

Below we review the main features of the integrability of sinh-Gordon 
model which will be important for us later. 

There exists a family of flat connections parametrized by the spectral 
parameter ^. We will also use the parametrization 

e = (5.8) 

The flat connection J = dw + Jw dw has the following components 

" ~ 2 V-^e-^T -^7 J ~ 2 V-^e^ ^7 J ^ ' 

Flatness is equivalent to the compatibility of the associated linear system 

(9^ + J^^ = a^ + J^^ = (5.10) 



which in turn is equivalent to the equation of motion (5.7). Another useful 
decomposition of the flat connection is 

J =^^y,dw + A + ^^yjdw (5.11) 
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using which we may write the string action as 

dxdx + dzdz 



2tr$^<l'^ 



(5.12) 



Certain specific solutions of the hnear system (5.10) associated with each 
puncture will be of key importance in the following. Since close to the punc- 
tures 

AV4 



and 
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1 



\ogT{w)T{w) 



we get two solutions, which close to the puncture behave like 



1 1 

_ 1 1 

W iWt 



'^9 = W ^^W 4^ 



1 1 

_ 1 1 

-W iVJi 



(5.13) 
(5.14) 

(5.15) 



It is clear that these solutions have nontrivial monodromies e^*^*^^^ around 
the puncture w = with 

r 



vr 



TT 



(5.16) 



It is in fact convenient to get rid of the vr by a gauge transformation = V'if 
with 



V 



^ I (w—wi){w—W2){w—w:j) 
{w—wi ) {w—W2 ) {w—w:j ) 







v 







{w—wi)(w—W2){w—w^) 
{w—wi)(w—W2)(w—W3) J j 



(5.17) 



Then our final basis of solutions associated to the puncture at w = wi is 



^1 
^1 



—j=\W — W\) 4« [W — W\) 4 ■ 

V2 



Ml 



u 



-1 



where the constants u\ are given by 



Ml 



lWl2Wi3j4 

(^12^13)4 
13 



(5.18) 
(5.19) 

(5.20) 



with Wij = Wi — Wj. The solutions 1 (i.e. and 1 (i.e. have the 
monodromies e*^^^-* and e"*^*^^^ with the pseudomomentum given by 

p(0 = A^(^e-^^ (^Avrsinh^^) (5.21) 

Several comments are in order here. These solutions can be continued to 
the neighborhoods of the other punctures. Of course we do not know their 
analytical expressions so we cannot perform this explicitly. Generically these 
solutions will no longer be eigenstates of the monodromy operator around 
other punctures. However, since the space of solution of the linear system is 
two-dimensional, we can express 1 and 1 as linear combination^ of an anal- 
ogous basis k and k at the puncture w = Wk- It is exactly these coefficients 
which are the key ingredients for the computation of the AdS part of the 
3-point correlation function. In order to fix an inherent ambiguity associ- 
ated with nontrivial monodromy, we have to fix once and for all the path of 
analytical continuation, whose detailed form will not be important for us. 

It is clear that the pseudomomentum of the elementary solutions obeys 
the important general property 

p(e-0 = -piO (5.22) 

This suggests that it should be possible to obtain the second solution 1 from 
the first 1. Since just changing C, — )■ e*'^^ modifies the expressions for the fiat 
connection, one has to perform in addition a similarity transformation 

UJ^,^{w, w; = J^d^, w- e^i) (5.23) 

with U = ia^ to compensate. Therefore the second solution can be obtained 
from the first ^'(w, ^) through 

^-k{w,w;0 = cr,^k{w,w;e'^0 (5-24) 

This is a crucial property which allows for the formulation of a set of func- 
tional equations for the overlap coefficients. 

Let us close this section with a comparision of the present set-up of a 
3-point correlation function in AdS2 with the case of Pohlmeyer reduction 
for polygonal Wilson loops in AdS^. 

"^with coefScients depending just on the spectral parameter 



14 



In both cases we have the same modified sinh-Gordon model, but with 
the modification defined in terms of functions with quite different analytical 
properties. In the case of polygonal Wilson loops we have a polynomial 
with a single asymptotic region (covered by several Stokes sectors), here we 
have a rational function with three (2"'^ order) poles and thus we have three 
distinct asymptotic regions. In the Wilson loop case, only the 'small' solution 
was unambigously defined, while here two solutions are uniquely specified 
as eigenfunctions of the monodromy operator at each puncture. Finally, the 
spacetime picture is quite different. In the Wilson loop case, the target-space 
was AdS^ and one had natural 'left-' and 'right-' linear problems. Here the 
target space is one dimension less {AdS2) and we have to develop appropriate 
reconstruction formulas and impose boundary conditions characteristic of 
a 3-point correlation function (i.e. fixing the boundary coordinates of the 
operator insertion points Xk and the target-space cut-off z = £). 



5.1 Overlaps 

In this section we will derive and solve functional equations for the overlaps 
between the elementary solutions associated with each puncture defined in 
the previous section. For any two solutions of the linear system (5.10) 
and one defines the antisymmetric product (skew-product) 

(kl) (5.25) 

which is the determinant of the matrix formed by the column vectors 
and It is a function of the spectral parameter ^ (or equivalently 6). Our 



elementary solutions (5.18)-(5.19) have the canonical normalization 

(kk) = 1 (5.26) 

A characteristic feature of the product {kl) is that for any four solutions 
the relevant products obey a purely algebraic relation called the Schouten 
identity: 

(tj) {kl) + {tl) {jk) + {ik) {Ij) = (5.27) 

In our case we have 6 distinguished solutions of the linear system - 
1, 1, 2, 2, 3, 3. Our aim is to find the skew-products between these solutions 
as functions of 6, given the set of conformal weights Ai, A2 and A3. 

It is convienient to repackage the products between the various solutions 
into connection matrices Mm which transform the coordinates of a solution 
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in the basis associated to the puncture / to the coordinates in the basis 
associated to the puncture k. 
The equation 




amounts to the following equality between the elementary solutions 



(5.28) 



This means that 



(5.29) 



(5.30) 
(5.31) 



Now taking appropriate products gives an expression for Mm in terms of our 
fundamental products. 



M, 



M 



-(kl) -{kl) 
{kl) (kl) 



(5.32) 



The obvious compatibility conditions between the connection matrices 

Mkm = MkiMim MkMk = id (5.33) 

are in fact equivalent to the full set of Schouten identities. This can be 
easily seen by considering various choices for the solutions entering Schouten 
identities and comparing with appropriate elements of the matrix products 
( |K33| . 

The full set of functional relations for the products (kl) thus comprises 
the compatibility conditions (5.33) and the vanishing of total monodromy 

fiiMi3fi3M32fi2M2i = id (5.34) 

where 









-iPk{S) 



(5.35) 



As they stand, the equations (5.33) and (5.34) are a comphcated set of 
constraints for 12 unknown products. The key property which allows us to 
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transform them into a set of solvable functional equations is the property 



(5.24). Using this construction we may relate 6 of the 12 unknown products 



to the other 6 but evaluated at a shifted value of the spectral parameter. 
Explicitly suppose that the k elementary solution (eigenvector of the mon- 
odromy matrix around the puncture Wk with the eigenvalue e*^*) is 



(5.36) 



Then the second solution k (with eigenvalue e *^'=) is obtained through 



(5.37) 



where the superscript '+' denotes the shift 9-^9 + i7r/2. Using the fac 



that a,. 



ttk and 6^ 



bk, it follows that 



[kl) 
[kl) 



(kl)^' 



(5.38) 
(5.39) 



Now the relations (5.33) and (5.34) become functional equations for just 6 
products. 



Solution of the functional relation 



In this section we will solve the full set of functional equations (5.33) and 



(5.34). Let us first define the three functions 



X32 


= (32) (32) 




= (32) (32) 


^21 


= (21) (21) 



(5.40) 

++ (5.41) 

(5.42) 

Once we determine them explicitly, the products (32), (32) and (21) will be 
expressed through convolution with a cosh kernel and zero-mode parts. The 
remaining products (21), (31) and (3l) will turn out to be expressed in terms 
of the first three and the pseudomomenta. In fact for our applications, it 
suffices to know the formulas for the products between the unbarred solutions: 



This follows from '^t^^^ = Xk'^k and an argument that Afe = 1. 
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(12), (23) and (31) - so all of them will be expressed through X32 and some 
permutation of indices. However in this section, for completeness, we will 
solve all equations. 

We start from the equation M32M21 = M31. This just expresses (31) and 
(3l) through the Schouten identities: 



(31) = 
(31) = 

Then define Yi and Is as 



(32) (2l)++ 
(32) (21)++ 



(21) (32) 
(21) (32) 



>3 



(32) (21) 

(31) 
(32) (21) 

(31) 



(5.43) 
(5.44) 



(5.45) 
(5.46) 



We can express (21) in terms of Yi. Plugging the results into the formula for 
Y3 we see that we can express X32 as 



32 



32 



Y, 



++I3 



(i-r++)(i-F3: 



(5.47) 



At this stage is convenient to rewrite the monodromy equation (5.34) in 
the form 

Ms2^2M2i = n^^Ms^n^^ (5.48) 



The entries of (5.48 ) are in fact the counterparts of the Y functions introduced 
by Maldacena and Zhiboedov ([52] fig. 5). 



Explicitly we obtain 



The equation (|5.48|) enables us to determine Yi and Y3 defined earlier 

I _ g«{P3-Pl-P2) 



Yi 
Y3 



g2ip2 g— 2«P2 

X _ g*(P3+Pl-P2) 

g2ip2 g— 2ip2 



(5.49) 
(5.50) 



Now we proceed to the equation M21M13 = M23 which enables us to 
determine X21, and finally the equation M13M32 = M12 determines X32. We 
check that all the remaining compatibility conditions of (5.33) are satisfied. 
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Consequently, the final functional equations for (32), (32) and (21) read 

P1+P2+P3 



(32) (32)++ 
(32) (32)++ 
(21) (21)++ 



sin P^-Pr^" sm ■ 



sinp2 sinp3 



sm 



P3-Pl-P2 



sinp2 sinp3 



sm 



P1-P2-P3 



sm 



P2-Pl-P3 



sm 



P2-Pl-P3 



(5.51) 
(5.52) 
(5.53) 



smpi smp2 

with the right hand sides being exactly our functions X32, X32 and X2i. In 
the above expressions we did not use the specific form ofpk{0) given by (5.21 ) 
but only the generic property 



Pk( 



ITT 



-vk{e) 



(5.54) 



Therefore, the above solution may have a much greater range of applicability 
than the specific case of no spin in AdS^ that we consider in the present 
paper. 

Let us now specialize to the pseudomomenta (5.21 ) and use the parametriza- 
tion ^ = e^. Then 

Pfe(6') = Afc7rsinh6' (5.55) 
The above functional equations can be recast in the form 

sinh( ^2+A3^^ cosh 6) sinh( ^i±^^±^7r cosh ( 



(32)^ (32)^ 
(32)+ (32)" 
(21)+ (21)- 



sinh( A27r cosh 6) sinh( A37r cosh 6) 
sinh(^i±^a^7r cosh 9) sinh(^i±^a^7r cosh ( 



45.56) 



sinh( A27r cosh 6) sinh( A37r cosh 6) 
sinh(^2±A3^7r cosh 9) sinh(^i±^a^7r cosh ( 
sinh(Ai7r cosh 9) sinh(A27r cosh 9) 



(5.57) 



(5.5J 



As mentioned before, for our purposes we will be interested in the solution 
of the first equation. The formulas for (12) and (31) can then be obtained 
simply by a permutation of the Aj's. 

The right hand side of the first equation has the property that it ap- 
proaches a constant when 9 — )• ±00 thus making the solution simpler. The 
basic functional equation to solve is 



f+f~ 

J a J a 



-aTT cosh i 



(5.59) 
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which can be solved by convolution 



fa{9) = exp 



de' log (1-6-'^^=°'^^^') 

2^ cosh(^ - 9') 



(5.60) 



Therefore we get the following expression for the product (32): 



(32) (6) 



le 



/a 2+A3-A1 (6')/a2+A3+Ai [d] 
f2AMf2AM 



(5.61) 



where the first term is a zero-mode part depending on two constants M 
and M*. These constants can be found from the leading WKB asymptotics 
of (32) (9) which can be found independently. We will discuss this part in 
section 8 and Appendix C. 



The formula (5.61 ) is the key formula of this section. We will use it in the 



following to obtain the AdS contribution to the 3-point correlation functions. 
Before we end this section, for completeness, let us discuss briefly the so- 



lution of equations (5.57) and (5.58). The right hand sides of these equations 



do not approach a constant when 9 — )■ ±00 so we cannot directly use the 



convolution with the cosh kernel. Apart from (5.59) we just have to consider 
in addition 



with the solution 



f+f- 

J a J a 



fa{9) 



62 



TT cosh I 



(5.62) 



(5.63) 



This will then solve the functional equations for (32) and (21). However in 
this case the zero mode part is undetermined. We will not consider this issue 
further since we do not need these expressions in the remaining part of the 
paper. 



5.2 Reconstruction formulas 

In this section we will show how one can reconstruct the string solution in 
the AdS2 target space from the Pohlmeyer data. The explicit expressions for 
the string solutions are important for two reasons. Firstly, the correlation 
function has to be regularized by making a cut-off ai z = £. This has to 
be translated into a worldsheet cut-off around each puncture. Secondly, we 
need to have control over the coordinates of the operator insertion points 
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Xk in gauge theory. In particular the standard conformal dependence on 
should arise automatically. 

We will show that the string solution can be reconstructed from two given 
solutions ^ A and ^'^ of the linear system for ^ = = 1) normalized by 
a^b) — 1- Equivalently, it is determined by the coefficients a, ft, 7 and 5 

of 

= a^i + *B = 7^'i + (5*i (5.64) 

satisfying aS — (3^ — 1. These two solutions can also be combined into a 
2x2 matrix as 

^ = ^) (5.65) 

Global embedding coordinates 

We will present reconstruction formulas in the global embedding coordinates 

Y^^^{l-x''-z^) Y^^^{l + x'' + z^) Y^^- (5.66) 
2z 2z z 



satisfying 

(y^)2-(y2)2+(y3)2 ^ _^ {dY^f-{dY''f+{dY^)'' = ""^ '/^ (5.67) 



dx"^ + dz"^ 
7i 



z^ 



Once we know y\ the Poincare coordinates may be easily extracted through 

y2 _ yi ^ 1 Y^ = - (5.68) 
z z 

Reconstruction formulas 

The reconstruction formula for the string solution is 

y^ = Itr (|a-^C^^D*) (5.69) 

where 

C={\'^ ^=(°n (6.70) 



^-1 oy \i 0^ 

and a are related to the standard Pauli matrices by 

a^^a^ a^^ia^ ^ (5.71) 
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Using the equations 



+ = 



+ = 



written in the original gauge (5.9), we may verify that 



(dY^BY^) 
ddY^ 



-1 

T{w) 
T{w) 



{dY^dY^)Y^ 



(5.72) 

(5.73) 
(5.74) 
(5.75) 

(5.76) 

(5.77) 



From the formula (5.69) we may now express the AdS2 coordinates di- 



rectly in terms of the components of \l/ given above: 
1 



= Y^-Y^ = 2iac 



X 

z 



Y^ = i{ad + bc) (5.7; 

Note that these expressions are invariant under the gauge transformation 

A 



A 



-1 



(5.79) 



thus we can use them also in our final gauge (5.17) 



The operator insertion points Xk and the target space cut-off S 

We may now use the above formulae to express the gauge theory operator 
insertion points in terms of the two solutions of the linear system \E'^, 
which determine the classical string embedding. Fortunately, close to 



the puncture we have explicit formulas (5.18)-(5.19) for the basis of solutions 
around each puncture. Using these formulas we see that for ^ = 1, the dom- 
inant solution around the puncture Wk is k. So only the (3 and S coefficients 



of \E'a,b in (5.64) will be relevant. 



Using the formulas (5.78) and the exphcit expression (5.19) we get the 



link between target space z coordinate and the worldsheet coordinate around 
the puncture w = Wk 



where /J^ = {k'^ a) 



(5.80) 



22 



This allows us to relate the target space cut-off ^ = £^ to the worldsheet 
cut-offs Ek'. 

Ak\ogek = \og£ + \og\{k^A)\^ (5.81) 

Similarly, we obtain expressions for the coordinates of the gauge theory op- 
erator insertion points 



(5.82) 



The two above expressions (5.81) and (5.82) are the key results of the 
present section which will be essential for the determination of the 'diver- 
gent' part of the AdS action integral in section 8 (recall also the overview in 
section 4 above). 



6 The AdS action 

After the above preparations we are now ready to tackle the calculation of 
the AdS contribution to the 3-point correlation function using Pohlmeyer 
reduction. 

We have to compute the action of the AdS2 part of the solution over the 
worldsheet, which is a 'regularized' 3-punctured sphere with 3 disks cut out 
around the puncture^ \w — Wi\ > Si. The worldsheet cut-off's around each 
puncture are not independent but are determined by the single target-space 
cut-off z = £ 

\f\ f dxdx + dzdz , , 

~v / — — ^^-^^ 

Let us emphasize that this is not the area of the worldsheet as there is a 
nonzero energy-momentum tensor. Using the elements of the Pohlmeyer flat 
connection the above integral can be written as 



— / 2tr$^<l>^ (6.2) 

Since in the above expression we have both an unknown integrand (i.e. which 
depends on the solution of the modified sinh-Gordon equation which we do 
not know explicitly) and an unknown integration domain (since the world- 
sheet cut-offs depend on the target-space solution), it is convenient, as out- 
lined in section 4, to split the integral into a cut-off independent finite piece 



^For the puncture at w = oo we define the worldsheet cut-off through |w| < I/Eo 
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which involves the unknown integrand but can be integrated over the whole 
punctured sphere and a cut-off dependent part with an explicitly known in- 
tegrand. 

— / (2ti^^^^- VrTcfw] + — [ VrTcfw (6.3) 

As outlined in section 4, we also adopt a similar regularization for the 5*^ 
part thus the cut-off dependent part will appear in the final answer with 
coefficient 2. We will evaluate the first integral in section 7, and the second 
integral in section 8. 



7 The regularized Pohlmeyer contribution 



In order to evaluate the first integral in ( |6.3[ ), we will proceed as for Polygonal 
Wilson loops, and pass to a gaug^ where the part of the fiat connection 
is diagonal 

<^^^W^^W-' (7.1) 

Fortunately it turns out that the diagonalized does not depend on the 
unknown Pohlmeyer function 7 and is expressed as 




$-=l ./rl (7-2) 

The diagonal components of become more complicated 



(T,_ _ 4v'T 4v'T _ - 2 V -1 '^uKiiT 2 V -1 «""i7 



4VT 4VT 

however, the important observation made in [5T| is that the diagonal com- 
ponents of each expression can be treated as a single function defined on a 
double cover S 

2/2 = T{w) (7.4) 

of the worldsheet S. 



^By this we mean redefining the solution of the linear system ^ = W'^ with some given 
matrix W depending on the worldsheet coordinates w. 



24 




Figure 1: Cycles on a genus 3 surface. For concrete computations it is 
convenient to make all the Bi cycles to pass through the same point on the 
last cut e.g. w = 0. 



In this manner one can rewrite the integral 



2tr<l>^<l>^ 



/TTd^w 

as an integral over S of a wedge product of two closed 1-forms 
/ (^2tr$^$^- VrTd^w 

with 



I 
2 



UJ /\r] 



UJ 



VTMd- 



1 



w 



2V^H(cosh7-l)rf^ + J-^ 



(7.5) 



(7.6) 



(^7) dw 
(7.7) 



where for simplicity we used the original Pohlmeyer function (see (5.5)). The 
dw component of rj does not influence the integral but is chosen so that r] is 
also closed {drj = 0). 

If S had genus g (which is the generic case for Polygonal Wilson loops), 
one would use Riemann bilinear identity (or reciprocity) to reduce the integral 
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Figure 2: The polygon is the standard representation of a genus 3 surface 
whose boundary is the curve Lg=^ composed of the cycles Ai and Bi each 
traversed twice. The infinitesimal circles C± surround the singularities of rj. 
Pq is the (arbitrary but fixed) base point for constructing the function F 
such that oj = dF. 

to products of integrals over cycles 




However in our case E has genus 0, and the 1-forms may have singularities at 
8 points (two copies of the 3 punctures and 2 branch points of the covering 
y'^ = T{w)). One possibility to proceed is to prove an analog of Riemann 
reciprocity directly for this case. The resulting expressions are, however, 
quite messy. In the end, we decided to adopt a slightly different strategy 
by treating the punctures as infinitesimal branch cuts and using Riemann 
reciprocity for a genus 3 Riemann surface (see figure [T| with an additional 
treatment of the singularities at the zeroes of T{w). Let us note that the 
1] 1-form is neither holomorphic or antiholomorphic and generic textbook 
formulas are not directly applicable. 

The idea of the derivation of the Riemann reciprocity formula is to rewrite 
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one of the forms as an exact form: 



uj = dF 



(7.9) 



where F{P) = jp^^ u. This can always be done on the Riemann surface minus 
some contour. Then one transforms the surface integral into a 1-dimensional 
integral over that contour using Stokes theorem: 



00 A r] = I uj A 7] = I d{Fr]) = / Ft] 
s Js\L Jt,\l J l 



(7.10) 



In our case, since uj is regulaijjat the zeroes of T{w), the contour L may be 
taken to be the sum of the standard contour for a genus 3 surface Lg^^ and 
two infinitesimal circles C± around the zeros of T{w) aX w = ±ia as shown 
on figure [2j 
Therefore 



UJ Arj 



Fr]+ Fr]+ Ft] 

Lg^S J C+ J C- 



(7.11) 



The first integral gives directly the standard bilinear expression 



3 



Lg=3 i=i JAi JBi J A, Jb, 



(7.12) 



where we used the fact that 



u = 



(7.13) 



'C+ JC- 

Let us now concentrate on the remaining two terms. Now, 



Ft] 



c+ 



C+ \JPo J JC+ l\JPo Jp+/ 

OJ f 7]+ f ( [ uj^ri= [ CO [ ri-inl (7.14) 

Po Jc+ Jc+ \Jp+ / JPo Jc+ 



6 



Recall that we are on E. 
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where the last integral is computed in Appendix B.l. Then, adding a similar 
expression for the second zero we get 

p+ f f 1 

CO V + ^ 1] - 2- iix- 

Po Jc+ JPo Jc- " 

uj I ri-2-nx- 

p_ Jc+ 6 

-2-i7r- (7.15) 
6 



+ / ]Fv 

c+ Jc- 



where we used the fact that the 1-form r] is regular everywhere apart from 
the zeros oiTiw) and so 



+ / 77 = 

c+ Jc^ 



(7.16) 



Moreover, one can even show that r] = (see Appendix B.2). In this way 
we arrived at the final equality. 



Further, inserting this result into the integral (7.11), and computing the 



and explicitly computing the integrals of the 1-form u 



periods in ( |7.12[ ) using again the regularity of rj outside the zeros of T{w) 

-27ri = , / UJ = 27ri=^ (7.17) 



A2 



A3 



we find that 



u /\rj 



E / ^ , 

i=l J A, JB 

A 

T 



f] — 2 ■ in 



6 



2m 



Bi 



B2 



V + 



A. 



V 



Bz 



6 



(7.18) 
(7.19) 



The integrals over the cycles Bi may be expressed by integrals between the 
punctures. From Fig. [T] and the antisymmetry of rj under changing of the 
Riemann sheet we find 



2 / r/ 

'C_ii 



7]+ I r], 

Jb2 



V 



B2 



2 I r, 

'Cla 



r/- / T] 

Bi JB3 



(7.20) 



(7.21) 
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Hence, 



'T, 



TT TT 

6 ~ 2 



({A^ -2A) f 7]- 2A^ ! ^ (7.22) 



u Arj = —2i 

Therefore the regularized Pohlmeyer contribution becomes 

j (^2tT<i>^<i>^- ^JT{w)T{w)d\^ = ^-^ (^(Aoo-2A) j r] - 2A^ j r]^ 

(7.23f 

At this stage we have reduced the computation of the regularized Pohlmeyer 
contribution to the evaluation of the integrals of the 1-form t] between the 
punctures. This cannot be done directly, as we do not know of course the 
explicit form of the Pohlmeyer solution 7 or 7. However, as shown in |5T] . 
the integrals of 1] can be related to the 6 — t- —00 [C, — )■ 0) asymptotics of the 
parallel transport of a solution along the curve which is a WKB line [53]. 

The main idea is to apply the well-know semiclassical methods where 
with the role of the Plank constant is played by the spectral parameter ^. 
Then, the linear problem 

{d + J)^ = (7.24) 

can be approximately solved with the leading contribution coming from the 
part 

^ ~ e'^^-^'^^^'^'^ (7.25) 

Clearly, the approximation is the best once we are on the WKB line defined 

as 

'1 



Im \^-^/T{w)wj = (7.26) 

For our purposes, however, it is crucial to know also the subleading term 
related to the part of the fiat connection 



(7.27) 



where fj is the 1-form r] without the subtraction term \/T{w)/2, exactly as 
it shows up in ^^i. 

The basic object we want to compute in this limit is the skew product between 
two solutions at the punctures j, k. Then, the prescription is the following: 
(i) take the known solution \l/j(w^) at w = w'j in the vicinity of the puncture 
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(ii) transport this solution via the parallel transport equation along a curve 
given by the WBK line equation to a point w = w'l^ near the puncture at Wk 
taking into account the leading as well as the subleading terms 

(iii) compare the resulting with the known solution ^'^(w^) at w = w' 
The resulting formula reads 



k- 



lim (jk) 



5 \/T{w)dw+^ log(«,,-t«^)+% log{u,fc-u,;^,) 

3 

\ Q y/7liS)dw+^ log(iDj-tD;) + % log(tDfc-tD^) 



(7.28) 



■ e^'^j ^(7.29) 



where the logarithmic terms are due to the exactly known form of the solu- 
tion near the punctures. Moreover, these subtractions render the expression 
finite and therefore allow to extend the integration exactly to the punctures. 
This formula may be now compared with the exact expression for the skew 



product at any ^ (5.61), which contains two undetermined zero mode con- 
stants M, M* . Fortunately, they are given by the first two terms of the WKB 
approximation. Then, the path integral of r; may be given by a combination 
of the —7- — oo asymptotic of faid) function. 

In this way we obtain the following explicit expressions for the period 
integrals 



7] 



C- 



/i(2A - Aoo) + /i(2A + Aoo) - 2/i(2A) 
/i(Aoo) + /i(2A + Aoo) - /i(2A) - /i(2Ac 



where 



h{a) 



/ — cosh t' log (1 — e 

J — oo ^ 



■a-K cosh I 



) 



(7.30) 



(7.31) 



Together with (7.23) it gives our final explicit expression for the regularized 



Pohlmeyer contribution. 



7.1 Comparision with numerics 



Since the above derivation of ( 7.23 )-( 7.30) was quite comphcated and in- 
volved many new ingredients, we decided to test the result by numerically 



30 




A 




numerics 


our formula 


0.2 


0.3 


0.04536 


0.0450779 


0.5 


0.9 


0.107649 


0.107622 


1. 


1. 


0.426311 


0.426166 


1. 


1.05 


0.429572 


0.429503 


2. 


2. 


0.517689 


0.517688 


2. 


3. 


0.488985 


0.488985 


4. 


4. 


0.523584 


0.523584 


4. 


7.99 


0.0152435 


0.0152435 



Figure 3: Regularized action density for A = Aqo = 4.0 and a compari- 
sion between the numerically evaluated regularized action and the analytical 



results of the formulas (7.23)- (7.30). 



solving the modified sinh-Gordon equation on the 3-punctured sphere and 
directly computing the regularized Pohlmeyer integral from the numerical 
solution. We give some details on the numerical setup in Appendix A, while 
here we just summarize the results and the comparision with the analytical 



predictions ( 7.23 )-( 7.30). 



In figure [3] we show a plot of the integrand entering the regularized 
Pohlmeyer action 

- (e2^(-'^) + T{w)T(w)e-^^^'"'''^) - ^JT{w)T{w)\ (fw (7.32) 



in a qudrant of the (compactified) complex plane. The upper and right 
borders are mapped to the puncture ai w = oo, while the puncture a.tw = +l 
is right in the middle of the lower border. The solution in the remaining three 
quadrants follows by symmetry. In the table we have shown a comparision 



of the numerical evaluation of (7.32) together with the analytical results 



following from (7. 23)- (7. 30). The numerics becomes more difficult and less 
reliable for small A's. In particular the deviations in the first rows of the 
table are within numerical errors estimated by changing the number of points 
of the numerical grid. The remaining results show excellent agreement with 
the analytical formulas following from periods and solutions of the functional 
equations for the overlaps. 
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8 The regularized divergent contribution 

In this section we will deal with the remaining contribution to the action - 
an integral of v^TT on the punctured sphere with specific cut-off's around 
each puncture \w — Wi\ > Si. Since we adopt a similar regularization for the 



S contribution (see (4.4)), we have in fact two such contributions 



exp<(-— / \/T{w)T{w)d'w} (8.1) 



The above integral can be evaluated explicitly in the small limit relevant 
for us. To do it, we note that it can be expressed as an integral of the wedge 
product of two closed 1-forms: 

/ Jt{w)T{w) (fw = const ■ / Vfdw A Vfdw (8.2) 

and transformed into a product of residues of \/T times appropriate (regu- 
larized) integrals over intervals between punctures. The outcome is 

/ Jt{w)T{w) <fw = Finite - ^A^log^oo -^AHoge- '^AHoge 

J^\{e,} 2 2 Z 

(8.3) 

where 

Finite = ^A^ (^2 log 2 + (1 + a^) log 2 + (-2 - + 2^1 + a^) log + 

+ (1 + a^) log(l + a^) - 4Vl + a2 log(l + Vl + a^)^ (8.4) 

Let us first concentrate on the logarithmically divergent part. Taking into 
account the coupling-constant dependent prefactor, and generalizing slightly 
to three generic anomalous dimensions, we get 

exp|v^(A^logei + A^loges + A^log^s)} (8.5) 



We may use the relations (5.81) to express the worldsheet cut-off in terms 



of the physical target-space cut-off z = £ and the products between the 
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elementary solutions 1, 2, 3 and one {"^a) of the two solutions appearing in 
the reconstruction formulas of section 5. We get 

exp {J2 log ^ + ^1 log I (1^) I' + ^2 log I {2A) P + As log I {3A) 

(8.6) 

where Aj = V^Aj is the unrescaled anomalous dimension. 

We will now express the scalar products (kA) in terms of the gauge theory 
operator insertion points Xk and products between the elementary solutions. 
The two solutions of the linear system \E^^ and which determine the 
target-space string embedding are completely specified by their coordinates 
in e.g. the 1, 1 basis: 

= al + (31 = 7I + SI (8.7) 

where aS — P'y = 1. Similarly we can express the 2 and 3 elementary solutions 



entering (8.6) in terms of 1 and 1: 

2 = kl + ll 3 = ml + nl (8.8) 

where k, I, m and n are appropriate overlaps evaluated at = (^ = 1). In 



terms of the above quantities, the part of (8.6) depending on the products 
becomes 

J2 log I {kA) P = Ai log A2 log (fc/5-/«)'+ A3 log {mf3-na)^ (8.9) 



Now we may use formula (5.82) to relate the hitherto unknown coefficients 



a, (3 and 7 to the operator insertion points: 

_ {3^b) _ m6-n-f 

Solving these equations with the constraint a6 — (3'-/ = 1 yields 

[3' = ^ (8.13) 

Im — kn X12X13 

{kl3-laf = -ilm-kn)-^^ (8.14) 
n X12X23 

(m(3 — naV = —ilm — kn)- — (8.15) 
I X1ZX23 
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In the above formula / = (12)q = (12)g^Q, n = (13)q, while 

lm-kn= (12)o (3l)o - {21)^ {13), = (32), (8.16) 
using Schouten's identity. 



Plugging the above expressions into (8.6), we obtain finally the standard 



CFT spacetime dependence of the 3-point function 

1 



^£12 j^l+^2-A3 ^^-jAi+A3-A2 ^^^A2+A3-Ai 



(8.17) 



multiplied by an additional contribution coming from the products of the 
elementary solutions 

g-5^pVA((Ai+A2- As) log(12),) + (Ai+A3-A2) log(13>o+(A2+A3- Ai) log(32)(,) -j^g-j 

Going back to our solution of the functional equations (and returning to the 
symmetric case of Ai = A2 = A and A3 = Aqo) we see that the products 
have the following structure at ^ = 0: 

(12)q = e^-"+^-ii ■ e^-" (8.19) 



where 



with 



(13)o = (23)o = e'^'i-+*^^- ■ e^^- (8.20) 



K^n = A;(2A - Aoo) + A;(2A + Aoo) - 2A;(2A) (8.21) 
i^ioo = A;(Aoo) + A;(2A + Aoo)-A;(2A)-A;(2Aoo) (8.22) 



, rf^log 1 -6-°'^^°^'^^) 

k(a) = / ^ — ^ 8.23 

^ ^ ' 27r cosh e ^ ' 



The zero-mode constants M_ii and Mioo are evaluated in Appendix C and 
after substituting into (8.18), it turns out that they exactly cancel the finite 
term (8.4). Thus the remaining contribution to the OPE coefficient becomes 
finally 

exp I [(2A - Aoo)iC-ii + 2Aooi^ioo] } (8.24) 
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9 The final AdS contribution to the OPE co- 
efficients 

We may now sum together the two contributions to the OPE coefficients 



(7.23)-(7.30) coming from the regularized Pohlmeyer integral and (8.24) 



coming from the regularized divergent integral. Both contributions have the 
same structure yielding 



riOPE 
^AdS 



exp 



(2A - Aoo)P-n + 2AooAc 



(9.1) 



with 



P_ii = /i(2A - Aoo) + /i(2A + Aoo) - 2/i(2A) 
Pioo = /i(Aoo) + /i(2A + Aoo) - /i(2A) - /i(2A, 



(9.2) 
(9.3) 



where 



h{a) = -h{a) — k{a) 



1 

2^ 



sinh 6 
cosh^^ 



log (1 - 



-air cosh ( 



) d9 (9.4) 



Several comments are in order here. Firstly, the above expression does 
not depend on any details of the operators entering the OPE coefficient 
apart from their anomalous dimensions, thus it is universal for this class 
of operators. Secondly, the above expression has to be supplanted by the 
regularized contribution 



y/X r 

= e 



contribution- 



'TT 



(9.5) 



so one cannot draw conlusions on the behaviour of real OPE coefficients 
qOpe _ (jOPE _(jOPE ^ since the latter part is currently unknown. Thirdly, the 

factor exp(— V^/6) seems quite surprising, however its presence is essential 
for sensible extremal and small Aj limits which we will examine shortly. 

In this paper we have mostly considered the symmetric case of two equal 
anomalous dimensions. It should not be difficult to extend these considera- 
tions to the generic case of three distinct anomalous dimensions. Repeating 
e.g. the analysis of the regularized Pohlmeyer contribution suggests the fol- 
lowing structure. 
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Let us introduce the parameters a,: 
ai = A2 + A3-Ai a2 = Ai + A3-A2 as = Ai + A2- A3 (9.6) 

Then the general answer should be 

C^;/ = exp |- + F{a,, a^, a,)^ | (9.7) 

where 

F(ai, ^2, as) = aih{ai) + a2h{a2) + ash{as) + (ai + 02 + as)h{ai + 0^2 + "3) 

- («! + a2)h{ai + 02) - («i + a3)^(«i + "s) - ("s + 0L2)h{as + ^2) (9.8) 

The structure of F(ai,a2,Q;3) is very similar to the structure of formula 
(7.11) in but with a different function h instead of a logarithm. Below 
we will see that (7.11) arises from our formula in the limit of small anomalous 
dimensions. 

Extremal limit 

In the extremal limit Aqo = 2A, all the terms with h{a)'s with nonzero 
arguments will cancel between each other leaving the term (2A — Aoo)/i(2A — 
Aqo). Now h{a) ~ — ^ (see Appendix D), so the remaining term will cancel 
with the —VX/Q giving the expected result 



C^;/(A,A,Ao,=2A) = l (9.9) 



Small A, limit 



For small arguments, h{a) behaves like (see Appendix D) 

h{a)r^-^-lloga (9.10) 
6a 2 

It turns out that contributions coming from the leading term will cancel out 
completely. The subleading logarithmic terms yield an expression 



.OPE,. . . ^ ^ ^(2A-Aoo)2A-A.,(2A + Aoo)^^+^-A^^-y 

(9.11) 



C^i;f(A,A,Aoo)^ 



(2A)4^(2Aoo)2A 

which coincides with formula (7.11) in 
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Large Aj limit and the Painleve transcendental 

For large arguments h{a) oc a* ■ e~'^", and thus their contribution is exponen- 
tially suppressed yielding a surprisingly simple universal limit independent 
of the conformal dimensions of operators: 



exp 



6 



(9.12) 



The simplicity of this result suggests that there should exist a much simpler 
direct derivation of the above resultB This turns out indeed to be the case. 
In order to study the large A limit it is most convenient to study the 



modified sinh-Gordon equation in its original formulation (5.2) 



where 



dd'y 



T{w) 



TT sinh 7 



A2 

00 



4 (1 



w 



2\2 



(9.13) 



(9.14) 



The advantage of using (9.13) is that 7 — )■ around the punctures. Hence it 
would seem naively that 7 = would be a possible solution of the equations 
of motion. This is not the case, however, as due to our genericity assumption 



on the nonvanishing of (5.1), 7 has to have logarithmic singularities 



7 ~ ± log liLi ± 2a| 



(9.15) 



at the zeros of T{w). Nevertheless, in the large A limit, when 



generically very large, in order to minimize the string action (5.1), we expect 



to have an almost vanishing solution with two narrow logarithmic spikes 
around the two zeros of T{w). 

Let us concentrate on the neighbourhood oi w = ia and introduce a new 
coordinate through w = u + ia. Then (9.13) takes the form 



c?c?7 



A2 

00 



2a 



4 (1 + a 
\^ — 

C2 



2\2 



MVMsinh7 



(9.16) 



''Thanks to Pedro Vieira for asking this interesting question. 
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Redefining coordinates again 



(9.17) 



yields the standard sinh-Gordon equation 

ddj = sinh7 (9.18) 

In the large A limit the problem becomes rotationally invariant and after 
introducing a new variable 7 = 2U and R = 2r = 2\v\ we obtain the equation 
for a Painleve III transcendent normalized as in |5i 



+ If/' = isinh2f/ (9.19) 
li 2 



The coefficient of the logarithmic singularity (9.15) becomes 

Ur^±\\ogR (9.20) 

Now we may use the results of [54J to evaluate directly the regularized 
Pohlmeyer action. We may rewrite the contribution around w = ia a.s 



VTT(cosh ^-l)d^w = J (cosh 2U - 1) ^ (9.21) 

Note the Svr which comes from an angular integral corresponding to a 27r 
angle in the original u and w coordinates. We then use the substitution 
from [21] 

^ cosh2f/ = ^^^RF^{R) (9.22) 



2 RdR 

to get 

Stt f d R 

~ ~2 



This integral can be evaluated exactly using the asymptotic properties of 
Fc{R) established in [51]. At large R, Fc{R) ~ —R/i up to exponentially 
small terms, while for small R, RFc{R) — > 1/18. Thus the above integral 
evaluates to 7r/12. Now taking into account two such contributions and the 
prefactor of the integral (9.21), we arrive directly at our universal large A 
limit: 

e-^lT2+T2j = (9.24) 
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10 Summary and Outlook 



In this paper we have computed the universal part of the OPE coefficients 
of three heavy operators with no Lorentz spins. This contribution comes 
from the AdS2 part of the string cr-model, and has to be supplanted with the 
contribution of the part in order to obtain the full OPE coefficient of the 
relevant operators. 

Wc employed the methods of Pohlmeyer reduction, which have been pre- 
viously applied with great success to the case of null polygonal Wilson loops. 
It is interesting to notice that different aspects of the strong coupling physics 
of A/" = 4 SYM (gluon scattering amplitudes, anomalous dimensions through 
2-point correlation functions, OPE coefficients) may be expressed by the same 
(modified) sinh-Gordon equation with all differences encoded in the analyt- 
ical structure of the modification functions. Despite the similarities, the 
differences are significant, especially in the analytical structure and target- 
space reconstruction (as we are dealing with AdS2 instead of AdSs) which 
make the generalization nontrivial. As a cross-check of our results we made 
a comparison of our formulas with direct numerical solution of the modified 
sinh-Gordon equation. 

Unfortunately, as the OPE coefficients of operators dual to semiclassical 
spinning strings have not been previously calculated, in general there are no 
independent results allowing for testing our final expression. Even in the 
case were we have some information, like specific BPS operators with large 
charges, our lack of knowledge of the contribution precludes a direct check. 

However, there are two limits in which we may cross-check our formula. 
First of all, in the extremal case we find that the AdS contribution to the 
OPE coefficient does not have any semiclassical piece, as expected in this 
case. Secondly, when the anomalous dimensions are small ('medium'-type 
operators) our result is exactly equivalent to the Klose-McLoughhn formula 
obtained by a classical extremalization procedure of three point geodesies 
related to the three operators. It is reasonable to expect that such a point- 
like string/geodesic approximation to three point correlators is acceptable 
for not too 'heavy' operators, which do not generate an extended surface. 

An additional consistency check is provided by the derivation of the cor- 
rect CFT space-time dependence of the three-point correlators which arises 
from the regularized divergent part. 

It would be interesting to perform a comparision with the case of heavy- 
heavy-light correlators. However it seems that such effects as leading order 
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backreaction to the classical 2-point solution and corrections to vertex op- 
erators would have to be included in the heavy-heavy-light calculations in 
order for the comparision to be made. 

There are several directions in which our work could be further developed. 
Obviously, as we consider entirely the AdS contribution to the OPE coeffi- 
cient, one has to perform an analogous analysis for the part in order to 
obtain the full OPE coefficient of the 'heavy' operators. We expect that the 

contribution will depend on the particularities of the states in questions 
and not only on their conformal dimensions. Moreover, all selection rules 
should show up from the contribution. 

The part of the problem seems to be significantly more sophisticated 
for at least two reasons. The first reason is just technical, and should not 
pose too much problems. Namely, the simplest string solution rotates in 

C which amounts to the reduction of the cr-model action to three 
dimensional case. Then, one should consider the corresponding Pohlmeyer 
reduction, again with a prescribed nonzero energy-momentum tensor. Since 
we deal with a higher dimensional target space, it leads to a more complicated 
(but still integrable) equations. A more serious problem is connected with 
the classical wavefunctions of the external states which must be taken into 
account as well. In particular it is not clear whether their contributions to 
the 3-point function would cancel out with their contributions to the 2-point 
functions when constructing a normalization independent OPE coefficient. 
Currently we lack an appropriate formulation with definite regularization 
prescription. 

Further generalization would involve the computation of the OPE coeffi- 
cients of operators carrying charges related to AdS^ momenta, but here again 
a consistent treatment of the wavefunctions would be needed. 

Finally, it would be very interesting to develop an analogous framework 
for higher point correlation functions and identify the general form of the key 
functional equations for the overlaps in this case. 

Acknowledgments: We would like to thank Pedro Vieira and Amit Sever 
for initial collaboration, and Volodya Kazakov, Kostya Zarembo, Kolya Gro- 
mov, Tristan McLoughlin and Arkady Tseytlin for interesting discussions. RJ 
was supported by Polish science funds as a research project N N202 105136 
(2009-2012). 
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A Details of the numerical computation 



In order to solve numerically the modified sinh-Gordon equation with the 
boundary conditions relevant for 3-point correlation functions, we have to 
redefine the Pohlmeyer function since both the original version 7 and the 
modified one 7 are not suitable for numerics. 7 has logarithmic singularities 
at w = ±ia, while 7 blows up at the punctures. A convenient choice which 
is finite everywhere is 

73 = ^7 + ^ log + - ^ log (kl' + «') (A.l) 

Another difficulty is to find a parametrization of the 3-punctured sphere 
which would allow us to use a rectangular grid of finite dimension. By sym- 
metry, one can restrict oneself just to the upper quadrant Rew > and 
Imw > 0. In this quadrant we used the following mapping oi w — x + iy: 

s = t = A.2 

x+1 y+1 ^ ^ 

and considered the unit square s G [0, 1] , t G [0, 1]. The boundary conditions 
on the vertical edge s = follow from reflection symmetry, while those on the 
horizontal edge t — follow from the condition that 7 vanishes on the real 
line. This in turn is a consequence of the symmetry properties of the target 
space solution under a reflection y — > —y which ensures that dyX — dyZ — 
at y = 0. 

The upper and right edges of the square get mapped to the puncture at 
w — 00, while w = 1 is in the middle of the lower edge. We use Chebyshev 
spectral interpolation in the standard way in the t coordinate and separately 
in two subintervals s G [0, |] and s G [|, 1] in order to allow for nonanalyticity 
at the puncture w = 1 (here (s,t) = (^,0)). 

We use a Python interface to the PETSc library (petsc4py) to solve the 
spectrally discretized PDE. In order to get convergence we had to use an 
automatic differentiation package (ADOL-C with Python bindings pyadolc) 
to compute the Jacobian and use the LU linear solver from PETSc instead 
of the default iterative one. 
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B Evaluation of the integrals 

In the vicinity of the positive zero w — ia the 1-forms in the leading order 
read 

/ 1 dw /T-, \ 

u DVw - la dw, rj (B.l) 

16D [w — lay 

where 



4(l + a2)2 

Now, it is convenient to introduce a new variable — w — ia. Since p e C+ 
i.e., is located on an infinitesimal circle around w = ia, we parametrize it as 
t = ee**^, where e is an infinitesimal parameter and (p the angular variable. 
Then, 

r uj = 2VDe'' [\^''f''idcj)' + o{e^) = l^/Oe^e^''^ - 1) + o{€^) (B.2) 
Jp+ Jo 3 



Now, we may insert this result into the contour integral 

c2tt 

12 



i ^'V - e-^^^)d</.o(6°) = =^ + 0(6°) (B.4) 



B.2 ^^ry = 

Here we prove that the integral of rj around a zero of T{w) vanishes. Without 
loosing generality we assume a — 1 which fixes the position of zero w± — ±i. 
The cut is chosen to join the zeros. Now, consider an integral around w^. 
The integral contour can be deformed in the following way 



V=\ + + / +/ +/ + 

C+ \J X, — ll ./— ll,0Ol ./0Ol,ll Jll,x/ \Jx, — l2 ./— l2,0O2 ./002,l2 Is 

(B.5) 
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where the subscript denotes the Riemann sheet of the doouble cover E and x 
a fixed point on the cut. Now consider a single pair of integrals of the form 

/ r,+ I V (B.6) 

>/— ll.OOl J —10,002 



Since r] is proportional to vT, it differs on the two sheets just by a relative 
sign. So the sum of the two integrals is zero. Repeating this argument we 
arrive at 

ri = (B.7) 



'c+ 

C The zero mode parts 

In general, the zero-mode constants are given by the following expression 

r^'^ 1 A A 

Mi2 = lim lim / -^/f^dw+-^\n{wl-w[) + ^\n{w2-W2) (C.l) 
w[^wi W2^W2 J^'^ 2 4 4 

where Wi, W2 are the end points for the WKB trajectory (poles of T{w)) with 
corresponding anomalous dimensions Ai, A2. In our case 

J ^/Ti^dw = ^ J ^^f^dw = 
yiT^ln ( ^ ■ °'+"'+^^^/'+"'' ') -2ln(w + Va^ + w^) 



4 



= # ■ In ( ■ i+'+^S"^^^ ) - ■ ln(w + Va^ + w^) (C.2) 

One may easily verify that the first term gives the expected asymptotics 
at w = 1 and w = —1 (with the correct ^ factor due to the fact that 
y/l + a"^ = ^-), while the second term is responsible for the asymptotics at 
w = 00. 

For the punctures Wi = —l,W2 = 1 we get a cancellation of singularities 



in (C.l) and find the following finite result 

^1 ^00, 2A + A^ A. 

J^-ii = — In—— — m — : in 

2 4A2 - A2 4 2A - Aoo 4 ^ ' 
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In order to calculate Moo,i it is convenient to perform a change of variables 
w ^ -. Then Moo,i = Mo,i, where 



4 1 



dz Ar 



A 



Mo,i = lim ^ , hm J ;^y^aMz^ + ^ H^i - ^i) + t ^^(^^ 



Therefore, 

Mo,i = lim2j^^j=olini4^z2= 



A 1 (1+^ g^+Hv'l^y"^ ' 



>^ln(-^l) + fln(^2-4: 



(C.4) 
(C.5) 



-¥-ln(^ + V^) 
The singular part Bit z[ ^ zi — is 

•In + y^a2 + ^)+^ln(-zi) = ^-In (^^) +^ In(-^i) = ^(In2+i7r) 



Ac 
4 



while the finite part at z[ ^ Zi — reads 



(C.6) 



lim,/ 



-4 -In 



^->-zi=0 4 



2^ i 



1 / 



The part at 2;^ ^ 2:1 = 1 is 



lim,/ 



i^->.Z2 = l 4 



•In 



4 



•in 



^2 J 



+ )+f ln(l-4) = 



lim,/ ^-In 



4 + 1 



a2-4-+Vl+a2./"2 + 4 
^9 V 



4V 



%-.ln(^^+^a2+-^j+^lnK-l)+^m (C.8) 



The singular parts ~ ln(z2 ~ 1) again cancel and wc arc left with 

A 1^ - ^ . In (1 + V^?TT) + f ZTT (C.9) 

Finally we get 

Mo,i = ^(In2 + i7r) - f hi + f In ^ - ^ -In (l + v^^^) 

(C.IO) 

which may be simphfied to (recall that M^^i — Mq^i) 
M^,, = ^(In 2 + m) + f In 2 - ^ • In (^^) - f In (?^) (C.ll) 
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D Small a asymptotics of h{a) 

We divide the integral into two parts 

1 r sinh^ e 



h(a) = —l r:q_in(i -e"""™^'^^) (D.l) 
^ ' In cosh^ ^ ' ^ ' 



1 

TT 



oo 

an cosh t 



/»00 /"OO 1 

/ coshein(l-e-""^°^^^) - / — ln(l-e 

Jo Jo cosh e 



(D.2) 



The first integral was calculated in [55] and has the following asymptotics 
for small a 



cosh^lnfl - e-'^"^"^^^) = -— + o(l) (D.3) 

6a 







In order to compute the second integral we notice that 

ln(l - e-'^'^™^^^^) = In(acosh^) + o(l) (D.4) 
Then, the integral is finite in this limit and reads 



oo 1 POO 1 

^ — avrcosh^N ' 



, ln(l - e-""™^^'') = / — — ln(acosh^) + o(l) (D.5) 
cosh ./o cosh 

lna + o(l) (D.6) 



TC 

2 



Inserting these results to the integral we get at the leading order 

h = -— --Ina (D.7) 
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